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Light Hadron Weak Matrix Elements 

Laurent Lellouch a 

a Centre de Physique Theorique, Case 907, CNRS Luminy, F-13288 Marseille Cedex 9, France 

I review this year's developments in the study of weak matrix elements of light hadrons on the lattice, with 
emphasis on K°— K° mixing and K — > tvty decays. 



1. Introduction 

Weak processes involving light hadrons in gen- 
eral, and kaons in particular, have contributed 
significantly to our understanding of fundamen- 
tal interactions over the years. Most recently, 
the measurement of a non- vanishing Re e'/e in 
K -> 7T7T decays by KTeV and NA48 f| has 
provided unambiguous evidence for direct CP vi- 
olation. 

In constraining the Standard Model (SM), the 
physics of kaons is complementary to that of 
B mesons. This is clearly visible in Figure [j], 
where the types of constraints imposed by weak 
processes involving either of these particles are 
displayed. Kaons provide an important con- 
straint on the summit of the unitarity triangle 
through the measurement of e, the parameter 
which quantifies indirect CP violation in K — > tttt 
decays. This constraint requires a description 
of non-perturbative effects in K°—K° mixing, 
parametrized by Br. Lattice results for this 
quantity are commonly used in unitarity trian- 
gle fits. There are two new results this year for 
Bk , obtained with domain- wall fermions, by the 
CP-PACS H and RBC f§ collaborations. 

Rare kaon decays and other FCNC kaon pro- 
cesses also provide good probes of physics be- 
yond the SM. Here too the lattice is contribut- 
ing. Donini et al. have computed SUSY, AS* = 2 
matrix elements 0] and the SPQcdR collabora- 
tion have performed the first lattice determina- 
tion of the electromagnetic operator matrix ele- 
ment, (Ti a \Q+\K ) , which contributes to enhance 
the CP violating component of the Kl — ► 7r°e + e~ 
amplitude in SUSY extensions of the SM Hfzj. 

The need for a non-perturbative technique such 
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Figure 1. K and B physics constraints on the uni- 
tarity triangle in a world without errors (from 

as lattice QCD, to study the processes discussed 
above is clear: the simple free-quark picture of 
weak interactions is severely modified by the non- 
perturbative effects of the strong interaction. One 
of the most extreme examples of this is, of course, 
the AI — 1/2 rule in K — > tttt decays: de- 
cays in which isospin changes by AI = 1/2 are 
greatly enhanced over those in which the change 
is AI = 3/2. Using the operator product expan- 
sion (OPE) to separate short and long distance 
contributions, one finds that most of the enhance- 
ment, in a QCD based explanation, must come 
from non-perturbative QCD effects in the matrix 
elements, ((7r7r)/ = o,2|WA5 , =i|^°), 01 the effective, 
AS = 1 weak hamiltonian, Has=i PiPl- 

The verification of this A J = 1/2 enhancement 
from first principles, as well as the calculation of 
e'/e and more generally the study of non-leptonic 
weak decays are some of the big challenges still 
facing lattice phenomenology. Though the prob- 
lems these processes pose are many, much has 
been learned over the years: 

1) The direct study of K — > tttt decays requires 
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one to consider four-point functions. These corre- 
lation functions are statistically much more noisy 
than the two- and three-point functions encoun- 
tered in most phenomenological studies under- 
taken on the lattice. To reduce the noise, one 
can consider transitions where all three particles 
are at rest or even study K — > tt (and K — > 0) 
transitions. In both cases, one uses chiral per- 
turbation theory (%PT) to relate the quantities 
computed on the lattice to the physical matrix 
elements. 

2) The renormalization of the dimension-six, 
four-quark operators which appear in the effective 
weak hamiltonian is difficult on the lattice: 

• There is mixing with other dimension-six op- 
erators. Some of this mixing is the same as in 
the continuum. But with fermion formulations 
which break chiral symmetry explicitly, such as 
Wilson fcrmions, there is additional mixing with 
wrong chirality operators. This problem, how- 
ever, is well understood and relatively well con- 
trolled. 

• There is also mixing with lower-dimension op- 
erators. This gives rise to power divergences, pro- 
portional to inverse powers of a, where a is the lat- 
tice spacing. The subtraction of these divergences 
is numerically very demanding. One must also 
be very careful. As pointed out some years ago 
by Martinelli, perturbation theory is likely to fail 
for such subtractions. Indeed, the coefficients of 
these divergences may have non-perturbative con- 
tributions of the form e~ / dg ^ 9 ' = clAqcd, 
which are formally smaller than any term in a per- 
turbative expansion. However, when enhanced by 
a power divergence, they can give a contribution 
which is of the same size as the physical quan- 
tity that is being calculated. A second concern 
is the definition of power subtracted operators 
in the presence of discretization errors. Again, 
a discretization error of, say, O (clAqcd), on a 
linearly-divergent quantity with mass-dimension 
one may combine with this divergence and give a 
finite contribution of O (Aqcd)- 

Many tools have been developed over the years 
to address the problem of renormalization (for de- 
tails, I refer you to Stefan Sint's contribution): 

• Lattice perturbation theory, performed in 



terms of a continuum-like, renormalized coupling 
constant can, in principle, be used to subtract 
mixings with same-dimension operators. 

• CPS symmetry fll^Jl2|| , which is a discrete 
symmetry of many four-quark weak operators, is 
a powerful tool for classifying the mixings these 
operators may have. A CPS transformation is a 
CP transformation, followed by a switching trans- 
formation which changes s into d quarks and vice 
versa. This symmetry is of course broken, but 
only softly by terms proportional to powers of 
(m s — rrid). There are also generalizations of CPS 
in which different pairs of quarks are swapped. I 
will generically refer to all of these transforma- 
tions as CPS. 

• Chiral symmetry, in the form of xPT or °f 
the sytematic exploitation of chiral Ward identi- 
ties can be used to guide, for instance, the sub- 
traction of power divergences in K — ► 7r ampli- 
tudes for chirally symmetric lattice-fermion for- 
mulations JIl[] or the construction of effective 
weak hamiltonians in the case where the lattice 
fermions break chiral symmetry explicitly jjllfu} • 

• Chiral symmetry is not sufficient to deter- 
mine the logarthmically divergent renormaliza- 
tions required to make certain quark bilinears 
or quadrilinears finite. For these subtractions, 
a non-perturbative renormalization (NPR) tech- 
nique was devised fl5|| . It is performed with 
quark states and makes use of a regularization- 
independent (RI) scheme. It requires gauge fixing 
and the existence of a window Aqcd <C /i <C 7r/a, 
where \i is the renormalization scale, /i 2 is usu- 
ally taken to be the momentum squared of the 
quarks. Aqcd <C ^ is necessary to guaran- 
tee that the renormalization constants obtained 
are independent of the states used to calculate 
them, fi <C 7r/a keeps discretization errors in 
check. This technique can also be used to com- 
pute the coefficient required to subtract the mix- 
ing of wrong-chirality operators in calculations in- 
volving four-quark operators. 

• Another NPR technique has been developed, 
which makes use of a finite- volume, Schrodinger 
functional (SF) scheme fl^j . It is associated with 
a non-perturbative renormalization group scaling 
and overcomes the problems encountered with the 
NPR techniques of |fl5|| . It is gauge invariant, 
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requires no window and is mass independent. It 
has not yet been applied to four-quark operators. 

• Finally, to circumvent the mixing problem 
completely, the authors of [0 have suggested cal- 
culating the hadronic matrix element of the T- 
product of weak currents, T { JY(x) J^(Q)} for 
a <C |x| <C I/Aqcd on the lattice. The results 
would then be fitted to the continuum OPE, giv- 
ing directly the matrix elements of the four-quark 
operators in the desired continuum scheme. Only 
the renormalization of bilinears would be neces- 
sary. The drawback is the need for very small 
lattice spacings which makes this approach costly 
numerically. 

While much is already known, contributions 
are still being made to the problem of renormal- 
ization and mixing. Two new suggestions were 
proposed this year for obtaining Bk with Wilson 
fcrmions but without the subtraction of wrong- 
chirality contributions usually required [ |l8| , |l9| |. 
The one-loop mixing of three- and four-quark op- 
erators with same-dimension operators has been 
analyzed in the domain- wall formalism pp| , pl[ . 
The renormalization of quark bilinears and four- 
quark operators, for domain- wall fermions, has 
been performed non-perturbatively in the RI 
scheme |22 ] . The same renormalization, for oper- 
ators which have no power-divergences, has been 
performed at one loop for Neuberger fermions 




3) A third problem is that present day lattices 
are only a few fermi long. On such lattices, it is 
not possible to separate the hadrons in a multi- 
particle state into asymptotic states. 

4) A fourth problem is that the lattice method 
"only" provides approximate, euclidean correla- 
tion functions. This forces us to confront the 
"Maiani- Testa theorem" [g5| which, in one of its 
guises, states that euclidean correlation functions 
which describe processes involving two or more 
final-state hadrons in the center of mass frame, 
yield amplitudes in which these particles are at 
rest. Thus, an euclidean correlation function cho- 
sen to describe K — ► tttt decays will yield an am- 
plitude for K(0) — * 7r(0)7r(0), which clearly has 
the wrong kinematics, unless ttik = 2m n . These 
statements will be made more precise below. 



Both problems 3) and 4) can be addressed with 
%PT. For the first, %PT enables estimates of lead- 
ing finite- volume corrections. It helps with the 
second problem in that it provides a means of ex- 
trapolating results to the correct physical point. 
In its quenched or partially quenched versions, it 
further gives a handle on quenching errors. There 
are new results in %PT this year which are di- 
rectly relevant for lattice studies of weak decays 
of kaons and which I briefly review in Section |[ 

While xPT seeks to correct finite-volume ef- 
fects and problems resulting from the "Maiani- 
Testa theorem" , a new approach to non-leptonic 
weak decays was proposed this year. It uses the 
fact that lattice volumes are finite to circumvent 
the "Maiani- Testa theorem" and yield directly 
K — > 7T7T matrix elements from euclidean corre- 
lation functions p6[ . Similar issues are currently 
being investigated by Lin et al. [ p7[|28| ] . 

Before closing this rather lenghty introduction, 
I would like to say a few words about domain- 
wall (and Neuberger) fermions, since they are be- 
ginning to play a prominent role in the study of 
light-hadron weak matrix elements. For details 
I refer you to Pavlos Vranas' talk p9j . These 
two recent formulations of lattice fermions have 
a significant advantage over the more traditional 
Wilson and Kogut-Susskind formulations. They 
have a full chiral-flavor symmetry at finite lattice 
spacing. There is, of course, a price to pay. For 
domain-wall fermions it is a fifth dimension with 
N5 — > 00 sites; for Neuberger fermions it is the 
inverse square root of a large matrix. The ques- 
tion then becomes: how small an N5 or poor an 
approximation to the inverse square root can one 
take and still have sufficient chiral symmetry? Let 
me concentrate on domain- wall fermions, which 
are currently more commonly used for weak- 
interaction phenomenology. The assessment is 
somewhat mixed. The mathematical convergence 
in 7V5, which is exponential asymptotically, can be 
rather slow for the couplings and volumes com- 
monly used, though improvements can be made 
|3(],|3l]]. However, it appears that in practice, for 
currently used simulation parameters, the chiral 
symmetry achieved may be sufficient for studying 
the physics of quarks whose masses are around 
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that of the strange |32|-|34|. The criterion used 
in these studies is the residual mass, m res , which 
essentially measures the amount the bare quark 
mass has to be shifted away from zero in order 
to have physically massless quarks. At (3 = 6.0 
with Ns = 16 and on a 16 3 x 32 lattice, the au- 
thors of H find that m res (MS, 2 GeV) ~ 4MeV. 
Nevertheless, it should be kept in mind that 
the amount of chiral symmetry required will de- 
pend on the quantity studied. Indeed, the ex- 
ponentially small chiral symmetry breaking can 
combine with power divergences induced by this 
breaking to give effects which may no longer be 
considered negligible. This problem is certainly 
relevant for matrix elements of four-quark oper- 
ators whose mixings with lower-dimensional op- 
erators are often made significantly worse when 
chiral symmetry is broken explicitly. This issue 
is currently being addressed by the RBC collab- 
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Given the additional cost and potential diffi- 
culties of chirally- improved fermions, one is en- 
titled to ask where, in weak matrix element cal- 
culations, is exact or nearly exact chiral symme- 
try absolutely necessary? This question is all the 
more justified that one appears to be able to ac- 
cess the desired physics with Wilson or Kogut- 
Susskind fermions in most cases. These calcula- 
tions, however, are sufficiently complex and diffi- 
cult that it is not clear, a priori, whether the over- 
head associated with chirally-improved fermions 
cannot be offset, at least partially, by their im- 
proved chiral behavior, by the fact that they are 
non-perturbatively O (a)-improved, etc. It is thus 
important that all approaches be pursued. 

The rest of this review is organized as follows. 
In Section |^, I review two new studies in xPT 
which are relevant for lattice calculations of weak 
decays of kaons. In Section ||, I discuss K°—K a 
mixing. A brief description of the AS — 1 effec- 
tive hamiltonian is given in Section |J. Calcula- 
tions of K — > 7T7T matrix elements from K — > tt 
and K — > amplitudes are reviewed in Sec- 
tion ^|, including two new ambitious studies us- 
ing domain- wall fermions. In Section ^, I discuss 
some of the issues surrounding the calculation of 
physical K — > ttit matrix elements from lattice 



K — > 7T7T amplitudes. A new approach to the cal- 
culation of non-leptonic weak decays is presented 
in Section ^. Section || contains my conclusions. 

The focus of this talk is on recent lattice de- 
velopments. I will unfortunately not have the 
space to cover the results obtained by other meth- 
ods. Please see |$6| and references therein for 
descriptions of some of the other possible non- 
perturbative approaches. 

2. Chiral perturbation theory results 

There are at least two new studies in xPT this 
year which are directly relevant for lattice calcu- 
lations of weak decays of kaons. 

In the first, results for the one-loop, O (p 2 ) cor- 
rections to the K — > 7r and K — > 7T7t matrix el- 
ements of the electroweak penguin operators Qj 
and Q s (see Eq. (|l^)), which belong to the (8,8) 
representation of SU(3)l x SU(3)r, were calcu- 
lated |37j . The authors find that the K — > tttt am- 
plitudes, including O (p 2 ) counterterms, can be 
obtained from a study of the m 2 K , m 2 and p^ -pk 
dependence of K — » 7r matrix elements. They 
further find that chiral loops can give rather sig- 
nificant (+(27 ± 27)%) corrections to the O (p°) 
results for ( (tttt) j=2 | Q 7, 8 1 ■ These results are 
obtained in regular xPT. It would be interesting 
to see how they are modified in quenched (q) or 
partially quenched (pq) xPT. 

In the second study, it is one-loop corrections 
to the matrix elements of the (8,1) and (27, 1) 
operators for K —> it transitions at rest with de- 
generate s and d quarks and for K — > with 
m s ^ rrid, and the relation of these matrix ele- 
ments to K — > 7T7T amplitudes, which are studied 
in pqxPT pq |. Not surprisingly, they find that 
K — > 7r and K — > amplitudes are not sufficient 
to determine all O (p 4 ) couplings required to ob- 
tain K —f tttt matrix elements at this order. Fur- 
thermore, they observe that the chiral logarithms 
are typically large and can depend strongly on 
N sea . By using their results to guide extrapola- 
tions of lattice results for these amplitudes to the 
chiral limit, one can hope to obtain reliable re- 
sults for the O {p 2 ) octet and twenty-seven-plet 
couplings, which are interesting quantities in their 
own right and for which a number of phenomeno- 
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logical estimates are available. 



The basic ingredient is the chiral rotation 18 



3. K°-K° mixing 

The most general AS = 2 effective hamiltonian 
for K°—K° mixing can be written in terms of the 
following operators: 

01 = [s a d a ] V -A[S(3dp] V _A 

02 = [Sad a }s-p[spdf3]s-P 

3 = [s a df3]s-p[sf3d a } S -p (1) 

04 = [s a d a }s-p[sf3df3]s + p 

05 = [Sad/3]s-p[spd a ]s+P , 

where a and (3 are color indices. In the SM, only 
0i contributes; in extensions, such as the MSSM, 
the other operators may also be required. 

These operators have positive and negative par- 
ity components, and only the former contribute 
to K°—K° mixing. For instance, 



01 = 



VV+AA 



20 



VA 



and Ovv+aa determines the SM contributions to 
this mixing. As is well known, the explicit chiral 
symmetry breaking present in the Wilson formu- 
lation of fermions implies that Ovv+aa will have 
finite mixings with wrong chirality operators: 



Ovv+aa(h) = Z V v+AA{a^,g 2 ) 

x \o vv+AA (a) + ^A i (g 2 )0 l 



(3) 



with i = VV-AA, SS + PP, SS-PP, TT. This 
mixing is particularly bad here because Ovv+aa 
is subdominant in the chiral expansion. 

3.1. Two proposals for getting around mix- 
ing with wrong chirality operators 

As shown by Bernard and collaborators many 
years ago, CPS symmetry protects the parity odd 
components of the operators in Eq. ([!]) from mix- 
ing with operators of wrong chirality Q. Thus, 



Ova(p) = ZvA(an,g 2 )OvA{a) 



(4) 



So the basic idea behind these proposals 
is to relate the matrix element of interest, 
(K°\Ovv+aa\K°), to a correlation function 
where the only four-quark operator is Ova- 



exp [ia75T 3 /2] 



(5) 



It is implemented in two different ways. 

In the first, one considers twisted-mass QCD 
& 

C = t/j (D w +m Q + iwf 5 T 3 ) ip (6) 

+s (D w + mo) s , 

where Dw is the usual Wilson Dirac operator, 
tp is the (it, d) doublet and mo, fio and rriQ are 
bare mass parameters. In the continuum, this la- 
grangian would be equivalent to the usual, three- 
flavor QCD lagrangian. On the lattice, however, 
because of the Wilson term, it is not. The authors 
then suggest adjusting the renormalized mass m 
and twisted mass /i in such a way that the an- 
gle of the rotation, a, in Eq. (§) is tt/2. Then, 
the physical Ovv+aa is —20va in the twisted 
theory, which means that only a multiplicative 
renormalization is required. 

The authors of jl9| consider a Ward identity 
associated with an infinitesimal version of the ro- 
tation of Eq. (|): 

Z P (P(ti)0 vv+AA (/x; 0)P(t 2 )) = Z VA {a^) 
x [2m Z P J2(n(x)P(ti)OvA(a; 0)P(t 2 )) 

X 

- Z s (S(t 1 )OvA{a;0)P(h)} (7) 
-Z s (P(t 1 )O VA (a;0)S(t 2 ))]+O(a) , 

with P = J2s d^ 5 s, S = ^2 S ds, II = d-y 5 d - uj 5 u 
and m = Z A (0\d^\7T°) /2(0\P^7r°) . 

Correlators containing the operator S are 
exponentially suppressed for t\, —t 2 — * oo, 
since they are dominated by scalar instead 
of pseudoscalar contributions. Thus, the 
correlator J2 x ( Il ( x ) P ( t l) VA(Q)P(t2)) yields 
(K°\0 V v+AA(f4\K a ) without mixing. Ex- 
ploratory results for Bk using this method were 
presented at this conference 0. 

Both methods are generalizable to other ma- 
trix elements. Of course, none of this is neces- 
sary with fermion formulations which have a chi- 
ral symmetry, such as domain-wall, Neuberger or 
Kogut-Susskind fermions. 
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3.2. Matrix elements for K°—K° mixing 
beyond the Standard Model 

Donini et al. have computed the matrix ele- 
ments of all operators of Eq. (|l]) between K° and 
K° states [|| . They work with a tree-level, O (a)- 
improvcd Shcikholeslami-Wohlert (SW) action at 
(3 = 6.0 and 6.2, in the quenched approxima- 
tion. Their strange and down quarks are degener- 
ate. They renormalize the matrix elements non- 
perturbatively in the RI scheme, match them onto 
other schemes at one loop and run them at two 
loops. 

To quantify and subtract residual artifacts that 
might remain in chiral behavior of the matrix el- 
ement of 0i, they fit its dependence on mass and 
momenta to the form: 



8 f 2 
3JK 



-f3m K +lPK -PR- 



,(8) 



where the parameters a and (3 are pure artifacts 
and the dots include higher-order terms in the 
chiral expansion to which they are not sensitive 
numerically. Q Br is then simply 7. 

In analogy with the definition of Br = 
(K \Oi\K )/y K m K , they suggest the following 
normalization for the matrix elements of the other 
operators (i = 2, • • • , 5): 



Bi(jj) 



(K°\Oi(n)\K° 

ff< m K* 



(9) 



instead of the usual, vacuum saturation approxi- 
mation (VSA) normalization, which is tradition- 
ally expressed in terms of the quark-mass com- 
bination [m s (p) + md(n)] 2 . The problem with 
giving the matrix elements in units of their VSA 
values is that they are then used, in phenomeno- 
logical applications, with uncorrelated values of 
the quark masses, thus compounding the uncer- 
tainty on the matrix elements with those on the 
poorly measured quark masses. This problem ob- 
viously disappears with the normalization of Eq. 

The authors do observe some lattice spacing 
dependence in their results, which never exceeds 
two statistical standard deviations. They choose 



Table 1 

Matrix elements of the AS = 2 operators ofEq. ^) 
between K° and K° states as obtained in \qj in the 
RI scheme. Also given are the corresponding renor- 



malization 


group invariant (RGI) matrix elements, 1 


defined in J 






(Oi) 


RI = 2 GeV) 


RGI 


(Oi) 


0.012(3) GeV 4 


0.017(4) GeV 4 


(02) 


-0.077(10) GeV 4 


-0.050(7) GeV 4 


(o 3 ) 


0.024(3) GeV 4 


0.001(7) GeV 4 


(04) 


0.142(12) GeV 4 


0.068(6) GeV 4 


<e>5> 


0.034(5) GeV 4 


0.038(5) GeV 4 



to average their j3 = 6.0 and 6.2 results and keep 
the largest statistical error. Their final results are 
summarized in Table [j]. It is worth noting that 
at 2 GeV in the MS-NDR scheme, the non-SM 
matrix elements are typically 2 to 12 times larger 
than the matrix element of 0\. The result (Oi) 
in Table [j] corresponds to the following value for 
B K : 



B^ DR (2GeV) 



0.68(21) 



(10) 



1 The fit they 
version of Eq. 



tually perform corresponds to a rescaled 



3.3. Bk with domain-wall fermions 

As already pointed out above, chiral symmetry 
facilitates the calculation of Bk- Domain- wall 
fermions are thus a good candidate for such a 
calculation. The first study of Bk with domain- 
wall fermions was actually performed some years 
ago p0[ . Two new quenched results were pre- 
sented at this conference, by the CP-PACS Q and 
the RBC jlj collaborations. Both calculations 
are performed with Shamir's variant of domain- 
wall fermions and degenerate s and d quarks with 
masses approximatively ranging from m s /A to 
m s . They differ, however, in the gauge action 
employed. CP-PACS use an RG-improved action 
at two values of (3 = 2.6 and 2.9, corresponding 
to an inverse lattice spacing of 1.81(4) GeV and 
2.87(7) GeV, respectively, as determined from the 
p-meson mass. RBC used a standard Wilson 
plaquette action, with j3 = 6.0, corresponding 
to an inverse lattice spacing of roughly 2 GeV. 
Another important difference in the calculations 
is that RBC renormalize the matrix element non- 
perturbatively ^2j, using the techniques of (lj|, 
while CP-PACS match their results onto the con- 
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tinuum perturbatively, at one loop [|20f . 

CP-PACS perform a rather extensive study 
of the dependence of their results on fifth- 
dimensional and spatial sizes and on cutoff. At 
(3 = 2.6, they work with the following four lat- 
tices: (16-24-32) 3 x 40 x 16 and 24 3 x 40 x 32. 
At (3 = 2.9, they have the two lattices (24-32) 3 x 
60 x 16. To investigate the chiral properties of 
their AS = 2 matrix element, they consider the 
ratio 



Br 



(K°\s~f 5 d\0)(0\d~f 5 s\K°) 



(11) 



which should vanish, by chiral symmetry, in the 
limit rriK — > 0. They study Bp in the limit of 
vanishing quark mass, obtained by linear extrap- 
olation from finite quark mass. They find that 
extending the fifth dimension by a factor of two, 
from 16 to 32 points, does not reduce the devia- 
tion of Bp from zero which they observe. Going 
to smaller lattice spacing does not reduce the ef- 
fect either whereas increasing the volume does. 
The effect is of order -10 to -20% in units of Bp 
at m K on their (32-24-16) 3 x 40 x 16_ lattices. | 
They conclude that the deviation of {K a \Oi \K°) 
from zero at vanishing quark mass must be a 
finite- volume effect. This is supported by the fact 
that the volume dependence of Bp at finite mass 
increases rapidly as this mass is reduced. At mx, 
the reduction in Bp in going from the largest to 
the smallest volume is small, less than 2%. At 
am q = 0.01, it is approximatively 15%. 

They then study the mass-dependence of Bk, 
and interpolate to the kaon mass using the fol- 
lowing %PT-inspired functional form 



B 



K 



B(\ — 3c x am q hi(am q ) + b x am q ) , (12) 



as shown in Figure || for their 24 3 x 40 x 16 lat- 
tice. The physical point is obtained at half the 
strange quark mass, as estimated from the ex- 
perimental value of mx/mp. The value of Bk 
extrapolated to the chiral limit appears to be ap- 
proximatively 30% smaller than the value at rriK, 
which would help reconcile lattice results with the 
chiral-limit result for Bk obtained recently in a 



It will be slightly larger for Bp extrapolated to m q 
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Figure 2. Dependence of B% DR (2 Ge V) on bare 
quark mass, from CP-PACS /§/. 

large- N c approximation to QCD jll). However, it 
is important to note that finite-volume effects of 
the kind described above could significantly dis- 
tort the chiral extrapolation of Bk- If so, an ex- 
trapolation to infinite volume would be necessary 
to determine Bk reliably in the chiral limit. 

Finally, they study the dependence of Bk at 
rriK on lattice-spacing, spatial volume and fifth 
dimensional size. They find that these depen- 
dences are small. They obtain the following result 
for B K 



B^ DR (2GeV) = 0.575(6) 



(13) 



by fitting, to a constant in lattice spacing, the 
results that they obtain from their runs on a 
24 3 x 40 x 16 lattice at (3 = 2.6 and on a 
32 3 x 60 x 16 lattice at j3 = 2.9. Errors in Eq. 
( |l2| ) are statistical only. 

From their fixed lattice spacing calculation on 
a 16 3 x 32 x 16 lattice, the RBC collaboration 
obtain the following result for Bk- 



r>NDR 



(2GeV) = 0.538(8) 



(14) 



This result is lower than the CP-PACS re- 
sult. This difference may be due to the fact 
that RBC rcnormalize their results perturba- 
tively whereas CP-PACS do so perturbatively. 
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Figure 3. B^ DH (2GeV) with domain-wall 
(filled symbols), Kogut-Susskmd (open symbols) 
and Wilson/Sheikholeslami-Wohlert (stick symbols) 
fermions. All results are quenched except for the 
Nf = 3 result of OSU and all errors are statistical. 
PT stands for perturbative matching, NP for non- 
perturbative matching a la J7^ / and WI, for match- 
ing using Ward identities. The references for the 
different results are: CP-PACS RBC Jg/, BS 
^Ojj (with NPR, \iep-lat/9909wk , JLQCD (KS) j$J, 
KGS Mi OSU andNj = 3 J±£, KSGP M, JLQCD 
(Wil) UKQCD gj/, APE BGS~g|/. 

The non-perturbative matching coefficient used 
by RBC is 7% smaller than the corresponding 
mean-field improved, one-loop coefficient com- 
puted in po|| . 

3.4. Bk summary 

I summarize, in Figure |3|, the results for 
B% DR (2GeV) obtained by different collabora- 
tions with different fermion formulations. The 
size of the error bars on the Wilson results is 
mostly a reflection of the subtractions which are 
necessary to restore correct chiral behavior. The 
fact that different actions give different results at 
fixed lattice spacing is a sign that for some of 
the formulations at least, there is still some way 
to the continuum limit. In this limit, however, 
all formulations agree within one and some stan- 
dard deviations. The domain- wall results exhibit 
rather good scaling. They are also systematically 



lower at fixed lattice spacing. If this trend sur- 
vives further scrutiny of the continuum extrap- 
olation and of systematic effects, we may have 
to slightly revise our canonical estimate of Bk- 
However, for the moment, because it involves the 
most extensive study of systematic errors, the KS 
result of the JLQCD collaboration should still 
be taken as the reference result. To this one must 
add uncertainties due to quenching and the fact 
that the down and strange quarks are degener- 
ate in the calculation. Sharpe ]49J , on the ba- 
sis of qxPT and the preliminary "Nf = 3" OSU 
results Q, suggests an enhancement factor of 
1.05 ± 0.15 to "unquench" quenched results for 
Bk- He advocates an additional factor 1.05±0.05 
to compensate for the effect of working with de- 
generate quark masses. Here, we choose to keep 
his estimate of errors, but not use his enhance- 
ment factors so as not to modify the central value 
on the basis of information which is, for the most 
part, not provided by lattice calculations. Thus, 
we quote 



r>NDR 
K 



(2GeV) = 0.628(42) (99) 
B* LO = 0.86(6) (14) , 



(15) 



B K LU is the two-loop RGI S-parameter 



where B% LO 

obtained from B^ DR (2GeV) with N f = 3 and 
a s (2 GeV) = 0.3. Of course, future studies should 
investigate systematically the effect of including 
light flavors of dynamical quarks. 

4. K ^ irir decays: general considerations 

With all massive modes, including charm, in- 
tegrated out, the AS = 1 effective hamiltonian is 
given by 3 



n AS =i = -j=V ud V* s y>, + ry t )Q t , (16) 



'AS / j V 
j=l 



where Zi and j/i are short-distance coefficients, 
r = -V t *V td /V: s V ud and 

Qi = (s a up)v-A {upd a )v-A 
Q 2 = (su) V -A {ud) V -A 

Q3,5 = {sd)v- A ^{qq)v T A 

q 



'For a review see, for instance, 



9 



!4,6 



{s a df3)v-A ^^{qpq a )v^A 



(17) 



Q&,10 — ^ {s a di3)v-A^^e q {qpq a )v±A ^ 

Q 

with q — u,d,s. It proves useful to split the 
AS = 1 operators into a sum of operators which 
transform under irreducible representations of the 

1/2 3/2 

isospin group: Q t = Q % ' + Q i ' . 

Integrating out charm is questionable since 
m c ~ 1.3 GeV. Above the charm threshold, one 
has the following effective hamiltonian: 

2 

(l-T)X>(Qi-Q?) 



G F 

W AS=1 = ~^ VudVl 



10 



(18) 



where Q\ 2 = Qi,2fu - * c], the sum over g now 
runs over u, d, s, c and yi and 2j are given Nf = 4 
values. 

The role of the lattice is to compute: 
(Qi)i = ((^)i\Qi\K) (19) 

for the two isospin channels 1 — 0,2 and i = 
1, • • • , 10, where Qi, 2 in Eq. ( [l9| ) are to be un- 
derstood as Qi.2 — Q1.2 m the case of an active 
charm. ^| 

5. K — ► 7T7T from X — > 7r and X — > 

The idea here is to compute three-point func- 
tions to obtain matrix elements of AS — 1 op- 
erators between K and 7r states, where the light 
quarks typically have masses around m s /2, and 
then use xPT to relate them to the correspond- 
ing, physical K — » -kit matrix elements. This last 
step is usually performed at lowest non-trivial or- 
der. In the case of AI —1/2 transitions, there is 
an unphysical contribution to K — > tt matrix el- 
ements of (8, 1) operators which does not appear 
in its physical K — > irir counterpart. This contri- 
bution can be subtracted by considering K — > 
transitions Ol. 



4 The conventions adopted here for yie normalization of 
isospin amplitudes etc. are those of [ 



5.1. Electroweak penguins with quenched 
Wilson fermions 

As part of their study of AS — 2 matrix 
elements, Donini et al. || have computed the 
AI = 3/2, electroweak-penguin matrix elements, 
(Q7.s)2, which give important contributions to 
Re e'/e. The parameters of the calcu latio n are 



the same as those presented in Section 3.2 



They use leading order xPT to get the K - 
matrix elements in the chiral limit through 



(Q 



Vim 2 fa 



7,8/2 



lim £?5,4 

n„— >0 



(20) 



which they obtain by linear extrapolation of their 
K° — K° results calculated for degenerate quark 
masses, m s — rrid — m q . Here, f w — 131 MeV. 
They find, in the chiral limit and at 2 GeV in the 
NDR-MS scheme: f\ 



(Q 7 (2GeV))? DR /(m%U) = 1.7(3) , 

8.1(8) . 



(Qs(2GeV))? DR /(mlU) 



(21) 



These values are obtained by one-loop match- 
ing from their non-perturbatively renormalizcd 
results in the RI scheme. This matching is rather 
poorly behaved as it induces a 40% change in 
((57)2 and a 27% change in (Qs)2- Since all other 
results for these matrix elements are, to date, nor- 
malized by their vacuum saturation value and this 
value is not provided, I will not attempt a com- 
parison with Eq. j2l|). 

5.2. Matrix elements of the AS — 1 opera- 
tors from quenched domain-wall QCD 

At this conference, the CP-PACS Q and 
RBC j|,[55| collaborations presented preliminary 
results from calculations of the matrix elements of 
all 10, AS* = 1 operators. They use domain- wall 
fermions and work in the quenched approxima- 
tion. They are also considering the case of ac- 
tive charm. Both collaborations work on lattices 
of size 16 3 x 32 x 16, with a domain-wall height 
M = 1.8 and cutoffs a -1 ~ 2 GeV. In addition to 
a Wilson plaquette action at {3 = 6.0 used by the 
two teams, CP-PACS also perform the calculation 



It should be remembered that the use of different nor- 
malizations for these matrix elements will yield different 
values for the matrix elements in units of GeV 3 . 
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with a RG-improved action at j3 — 2.6. f\ In these 
calculations, a pseudoscalar meson composed of 
degenerate quarks of bare mass am q ~ 0.02 would 
have a mass close to that of the physical kaon. 
While CP-PACS renormalize their matrix ele- 
ments at one loop J21J] , RBC perform this renor- 
malization non-perturbatively [p2f . 

AI = 3/2 transitions involve no "penguin" con- 
tractions and therefore no mixing with lower di- 
mensional operators. And because of the approx- 
imate chiral symmetry of domain- wall fermions, 
the structure of the renormalization is the same 
as in the continuum, up to corrections which are 
exponentially small in N$. 

AI =1/2 transitions can receive contributions 
from lower dimensional operators. To leading or- 
der in the chiral expansion, the matrix elements 
1/2 

(7r + |Q/ (a)\K + ), i 7^ 7, 8, are quadratically di- 
vergent due to mixing with the operator 



S ub = (m s + m d )sd - (m s - m d )s"/ 5 d 



(22) 



To subtract this divergence, one uses the fact that 
the coefficients of the parity even and odd com- 
ponents of O su b are fixed by chiral and CPS sym- 
metry dllj . Then one constructs a subtracted op- 
erator Qi — UiO su \, by imposing the constraint 



{Q\Qi-aiO sub \K Q ) =0 , 



(23) 



to determine a>i. Because CP-PACS works with 
degenerate quarks, they obtain a, from a deriva- 
tive of Eq. ( p3[ ) with respect to m s at the point 
m s = ma- RBC work with m s 7^ rrid and obtain 
Ui from the slope in m s — rrid- 

Note that this construction relies heavily on the 
fact that the coefficients of sd and sj^d in O su b 
are fixed. In the absence of chiral symmetry, as 
when using Wilson fermions, this is no longer the 
case and one finds that the matrix elements dis- 
cussed above are cubically divergent due to mix- 
ing with sd and that the coefficient of this di- 
vergence cannot be obtained by studying K — > 
matrix elements. 

Using leading order xPT, CP-PACS relate the 
K — » 7r matrix elements that they calculate to 



6 For RBC, the gauge ensemble is the same as the one used 
in their Bx calculation. 



the desired K — > 7T7t matrix elements, through: 



<Q i )/ = C / /.(m|-m2)^ 



(24) 



-\Qt 



O sub \K^ 



I M70 75 d 1 0) (0 1 S70 75 u I K A 



[GeV J 



fife 

7 2 



(25) 



M7o75d|0)(0|s7o75w|-fr H 



\GeV d 



with i = 1, 2, ■ • • , 6, 9, 10, j = 7,8 and where 
AI = 1/2(3/2) and Cj = y^iVS) for / = 0(2). 
The Zij/Z\ are the constants required to match 
the lattice results onto the NDR-MS scheme. 
aJ 2 = and aj" 1 is obtained as described around 
Eq. ( p3| ) . Trip is the mass of the pseudoscalar me- 
son in the simulation. [. . .] x indicates that the 
term in brackets is linearly extrapolated to the 
chiral limit. 

Their results for the plaquette and RG- 
improved actions are very similar. We will there- 
fore concentrate on the former, because they can 
be compared directly to RBC's results, obtained 
with the same parameters. 

Figure || shows the subtraction of the power di- 
vergence in the calculation of (Qe)o performed by 
CP-PACS, as a function of their degenerate quark 
mass. The cancellation is severe and leads to a 
value for (Q|" 6 ) [| which is very roughly 40 times 
smaller than the individual contributions. That 
a signal remains at all must be due to the strong 
statistical correlation between the two terms. 

While CP-PACS presented preliminary results 
for the matrix elements, (Qi(2GeV)) I IDR , with 
i = ].,•■• 10, and / = 0,2, it is important to 
remember that this is a very difficult calcula- 
tion and that premature phenomenological con- 
sequences should not be drawn. I have therefore 
chosen not to quote these results here. Let me 
instead comment on the problems they may be 
confronted with. 



7 By Qf , I mean the operator with power divergences 
subtracted, but that still requires a logarithmic renormal- 
ization. 
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Figure 4. Subtraction of the quadratic divergence 
which appears in the construction of {Qe)o given in 
Eq. (pE^), for different values of the degenerate quark 
mass, m q . 

One of the features of their results is that some 
of the a s corrections in the matching of (Qt ub )o 
onto the NDR-MS scheme are large, as large as 
50% in some instance. It will be interesting to see 
what RBC finds with its non-perturbative match- 
ing. 

Also, as mentioned above, the fact that the di- 
vergence is quadratic in the AI —1/2 channel is 
a consequence of chiral symmetry. Because chi- 
ral symmetry is only approximate with domain- 
wall fermions, we actually expect there to be a 
cubic divergence due to mixing with sd, as in 
the Wilson fcrmion case. However, instead of be- 
ing of order 1, the coefficient of this divergence 
will be exponentially small, leading to a term of 
the form e~ cN5 sd/a 3 , which is formally of order 
m res sd/a 2 . Such a term will not be subtracted 
by the condition of Eq. (|23|). It will appear as a 
non- vanishing intercept at m q — —m res in a plot 
of the matrix element versus quark mass of an 
operator such as Qe, which should vanish in the 
chiral limit. Such an intercept is seen by RBC, 
as shown in Figure |[ This intercept is still rea- 
sonably small on the scale of the unsubtracted 
( 7r+ \Q6( a )\K + ) matrix element of Figure ||, at 
the values of am q at which the calculation is per- 




0.02 0.03 
am. 



Figure 5. K + — > tv + matrix element of the bare op- 
erator Qq. The line corresponds to an uncorrelated 
linear fit to the data. 



formed. However, it is large when converted to 
the scale of (QQ Ub )a in Figure || and increases like 
1/mp when rap is decreased, as the definition 
of Eq. ( p4| ) implies. [] Of course, the estimate 
of the size of this unphysical contribution will 
depend sensitively on how the chiral extrapola- 
tion of (TT + \QQ(a)\K + ) is performed, a delicate 
question in quenched QCD and in a finite vol- 
ume. ^] But this should be taken as a warning 
that such effects may be important and should 
not be ignored. One way to subtract such an 
unphysical intercept is to obtain the chiral-limit 
matrix element from the slope in the dependence 
of (ir + \QQ Ub \K + ) on m q . This is the procedure 
advocated by RBC J3^] and I would favor it over 
the use of Eq. (|J). 

It should also be remembered that what is cal- 
culated are K — > 7r matrix elements where the 
K and the 7r are degenerate and at rest. These 
matrix elements, suitably normalized, are then 
extrapolated to the chiral limit and translated 
into K — > 7T7T matrix elements using leading order 



8 Here I am assuming that the intercept measured by RBC 
is similar to the one CP-PACS would measure: the simu- 
lation parameters are identical. 

9 Note that the size of the intercept is not inconsistent with 
it being of order — 2m,. es (gg)/(a 2 /^), which is the chiral 
value of m res {ir + \sd\K + ) /a 2 . 
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XPT. This procedure thus requires a good control 
over chiral extrapolations, which is not trivial in 
the quenched approximation and in a finite vol- 
ume. Furthermore, the leading order x?T re- 
lations between K — » it and K — > tttt neglect ef- 
fects which may be important, such as final-state 
interactions. 

Finally, results involving an active charm quark 
have not yet been presented. Given the fact that 
the charm is not that heavy, its contribution as 
an active quark may be important. 

6. K — > tttt from unphysical K — > tttt 

This approach is complementary to the K — ► 7r 
and X — > method discussed above. Here, four- 
point functions are used to compute the matrix 
elements for K — > tttt transitions at unphysical 
values of the mesons' masses and momenta, im- 
posed by computational and theoretical limita- 
tions. Then, low order xPT is used to extrapo- 
late the result to the physical point. While this 
method is numerically more demanding because 
of the four-point functions required, renormal- 
ization is simplified. In fact, this approach ap- 
pears to be the only one possible when consider- 
ing AI = 1/2 transitions in the absence of the 
GIM mechanism (as is the case when studying e') 
with fcrmions which explicitly break chiral sym- 
metry, such as Wilson fermions. 

The SPQcdR collaboration are undertaking a 
quenched study of the AI = 1/2 rule and of e' 
using non-perturbatively O (a)-improved Wilson 
fermions f?J and NPR matching a la All 
mesons are taken at rest. They work with de- 
generate quarks, i.e. m u = rrid = m s , but con- 
sider also the situation where the quarks in the 
kaon and in the pions have masses tuned so that 
rriK — 2m n . With pions at rest, this latter situa- 
tion yields K — > tttt at threshold. 

In the A J = 3/2 channel, renormalization is 
particularly simple as there is no mixing with 
lower-dimension operators and CPS symmetry 
guarantees that there is no mixing with wrong- 
chirality operators of dimension six [ fL2|] . They 
presented an exploratory study of the matrix el- 



ement (Q2 + Qi)i, which constitutes the denom- 
inator of the AI = 1/2 rule. This amplitude 
has already been studied in the quenched approx- 
imation with unimproved Wilson fermion, per- 
turbative matching and degenerate quarks f3l| ~- 
|53| , as reviewed last year Q. By account- 
ing for finite- volume effects and chiral logarithms 
p5| , and using modern techniques in perturba- 
tive renormalization, the authors of the most re- 
cent study |33| were able to reconcile, with exper- 
iment, the physical amplitude obtained from the 
chiral-limit lattice results. The pioneering stud- 
ies of [51 5^| had found a discrepancy by a factor 
of roughly two. Non-perturbative matching and 
O (a)-improvement, as well as the situation where 
rriK = 2m n , as considered by SPQcdR, should 
help further clarify the situation for these decays. 

SPQcdR also presented exploratory results 
from the first K — > tttt calculation of the am- 
plitudes ((57,8)2 @- Combining their results with 
those obtained from K — > 7r amplitudes should 
help obtain more reliable predictions for these 
matrix elements. 

In the AI = 1/2 channel, SPQcdR are con- 
sidering mixing with the lower dimension opera- 
tors 



Op = (m s - m d )sj 5 d 

Og = {m s - m^sa^G^d 



(26) 



where the factor of (m s — wid) is required by CPS 
symmetry. (The AS = 1 operators of Eq. ( [Tt| ) 
are all even under CPS.) Clearly, working with 
m u — rrid — m s eliminates this mixing com- 
pletely |12|. It is shown in jl7j that the choice 
rriK = 2m n with all mesons at rest also simpli- 
fies the renormalization because no momentum is 
transfered by the AS = 1 operator. SPQcdR is 
pursuing both these avenues . 

7. Non-leptonic weak decays from finite- 
volume correlation functions 

Determinations of K — > 7T7t amplitudes, both 
from K — * 7r (and K — ► 0) matrix elements and 
from unphysical K — > tttt amplitudes, rely on low 



10 This issue is 
collaboration It 



currently being investigated by the RBC 



1 Og is actually subleading in the chiral expansion pq]. 
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order xPT. They therefore neglect chiral correc- 
tions which may be important. Aside from this 
limitation, it would be satisfying to be able to 
calculate K — > -kit amplitudes in full, directly on 
the lattice. Can this be done? 

7.1. Euclidean correlation functions and 
the "Maiani- Testa no-go theorem" 

The lattice method provides approximate esti- 
mates of euclidean correlation functions. While 
the Osterwalder-Schrader theorem guarantees 
that these correlation functions can, in princi- 
ple, be analytically continued back to Minkowski 
space, such a procedure is unstable in practice, 
when applied to approximate data. This led Ma- 
iani and Testa to investigate what can be ex- 
tracted from euclidean correlation functions with- 
out analytic continuation |p5| . 

They considered a correlation function which, 
in the case of K — » mr decays, can be written 

Ck-ktt {ti,h,t 3 ;p, -p) = (27) 

(7r^ 1 )^_ p -(^)H A s=i(0)^o(- < 3)) , 

where 



7Tj?(t) = / d 3 xe- l P- S Tr(x,t) , 



(28) 



and similarly for Kp{t). Then they investigated 
the behavior of this correlation function as a func- 
tion of t2 > in the limit t\, t$ — > +oo. 

So as to disentangle euclidean from finite- 
volume effects, they chose to work in the limit 
of asymptotically large volumes. This has for 
consequence that their spectrum of two-pion fi- 
nal states is continuous which, in turn, means 
that when i 2 is taken to +oo, only the ground 
state contribution can be picked out. Thus, 
they found that the euclidean correlation func- 
tion of Eq. (|27| ) only gave them the matrix el- 
ement (Tnr\HAS=i\K) with all mesons at rest, 
which is certainly not the physical kinematics. 
The further found that Ify/fa corrections to the 
t% — > +oo result contained information about the 
7T7T scattering length. 

This direction has been pursued by Ciuchini et 
al. who argue that they can reconstruct the de- 
sired weak decay amplitudes under the assump- 
tion that final-state interactions are dominated 



by nearby resonances and that the couplings to 
these resonances is smooth in the external mo- 
menta p3(|| . Lin et al. are currently working on 
eliminating the need for these assumptions ]27|| . 

7.2. Two-pion states in finite volume 

In 0], a different route has been followed. In 
present day simulations, lattice have sides, L, of 
order a few fermi. This means that the spectrum 
of two-pion states on such lattices is discrete: 
the quantum of momentum is Ap = 2tt/L = 
1.2 GeV/L[fm]. Thus, this spectrum is far from 
continuous. This statement can, in fact, be made 
much more precise. 

Consider a box of volume L x L x L with peri- 
odic boundary conditions and sides L > 3 fm. R 
In such a box, it was shown by Luscher f37| that 
two-pion energies, in the Af ( "spin-0" ) sector and 
center of mass frame, and below the four-pion 
threshold, are given by 



Wt = 2Jml + k 2 e , 



1,2, 



(29) 



£ir - 5f>(ke) = 4>(q e 



2tt 



where <5g(fc) is the irir scattering phase shift in the 
spin-0 and isospin-/ channel and <p(q) is a known 
kinematical function. 

To get some idea of what these equations say, 
one can expand their solution in inverse powers 
of L and one finds (£ < 6) 



Wt = 2y/ml + £(2tt/L) 2 + O (l/i 3 ) 



(30) 



This is the energy of two free pions, each with 
momentum \/I(2ir/L), up to corrections which 
go like one over the volume and which are de- 
termined by the phase shift. 

In Figure [| the two-pion energies are plotted in 
units of the pion mass as a function of the length, 
L, of the box's sides. The dashed curves are the 
energies of two free pions. The solid curves are 
the energies of two interacting pions with I = J = 
0, as derived from Eq. (29), with 6® calculated 
at one-loop in x?T 58 59j and taken from [|iio||. 



12 This guarantees that single-particle states resemble their 
infinite-volume counterparts up to terms exponentially 
small in L. 
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Figure 6. Two-pion energy spectrum in QCD below 
the inelastic threshold and as a function of the length 
of the box's sides. The dashed curves correspond to 
free pions. The solid curves correspond to interact- 
ing pions in the I — J = channel. Also shown is 
the kaon mass, mic, and the intersection of the first 
excited state (1=1) with this mass (dot). 

Please note that x^T is only used here for the 
purpose of illustrating how Eq. (^9|) works. It is 
not in any way required by the method. 

The first thing to note is that the spectrum is 
far from continuous for any value of L one may 
hope to reach in numerical simulation. The sec- 
ond point is that the corrections to the spectrum 
brought about by tttt interactions are not that 
large. This is due to the fact that they appear 
suppressed by 1/L 3 and the volumes considered 
are not that small. A third point is that Eq. 
( f2g| ) can be turned around. Thus, a numerical 
study of the finite- volume two-pion spectrum can 
be converted into a determination of the phase 
shift 6&(k) H . 

7.3. K — ► tttt decays in finite volume 

Having established that the two-pion levels are 
discrete, one can then imagine isolating K decays 
to some of these excited two-pion levels, thus de- 
termining the matrix element 

M/ = (St™)iI\H* s =i\K) . (31) 

13 A possible approach to isolating excited two-particle lev- 
els is discussed in fell. 



In Eq. (|T]) , all finite volume states are normalized 
to one: (K\K) = ((irir)i£\(inr)i£) = 1. 

Furthermore, having obtained <5g(fc) from a nu- 
merical study of the finite- volume two-pion spec- 
trum, one can imagine tuning L in such a way 
that the n th tttt state has energy W n = tuk- 
This is illustrated by a dot in Figure [| for n = 
1. The corresponding transition matrix element, 
M^, conserves both energy and momentum and 
the kaons and pions have their physical masses. 
It is not, however, a finite- volume matrix clement 
that we are after. What we want is the matrix 
element which describes the decay of a kaon at 
rest into a state of two asymptotic pions with en- 
ergy m K , i.e. = {(tttt out)/; m K \TLAS=i\K) . 
The hard work enters in showing that these two 
matrix elements are related. 

In p6[ , it is shown that the infinite- volume am- 
plitude, Tj{, is related to the corresponding finite- 
volume matrix element, M*, through [^] 

' R {q<P , (q) + k5 I '(k)} k=kn (32) 



\Tl\ =87rL b 



Mi 



This result assumes that W n < Am^ and that 
the final two-pion state |(7r7r) jn) is not degenerate 
(i.e. n < 8). It is accurate up to exponentially 
small corrections in L. 

There are various ways to proceed to prove this 
relation pq ]. One approach is to study the nor- 
malization of tttt wavefunctions in finite and in- 
finite volume, with the quantum-mechanical for- 
malism of [^7]|32|]. Then one uses the fact that 
the transition matrix elements probe the S'-wave 
component of the two-pion wavefunction near the 
origin and that this component differs in finite 
and infinite volume only in its normalization and 
possibly phase. 

Instead, one can switch on the weak interaction 
(Was=i) an d compute its influence on tttt ener- 
gies. This can be done directly, using ordinary 

14 In the study of the normalization of tttt wavefunctions 
in finite and infinite volume, as well as in the perturbativc 
check, both of which are discussed below, a relation which 
hold s fo r W n 7^ mjc can also be obtained. It is identical to 
Eq. (p2|), except that the factor of (m^/fc n ) 3 is replaced 
by W^nifc /k n . Whether this result holds more generally 
has yet to be verified. 
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quantum-mechanical perturbation theory, or one 
may start from Eq. (^9j), taking into account the 
effect of the weak interaction on the scattering 
phase. A comparison of the two results yields the 
relation of Eq. (§§). 

One can also verify this relation by working 
out both finite and infinite-volume amplitudes in 
perturbation theory, in a low-energy effective the- 
ory such as xPT. Because Eq. (32) does not de- 
pend on the details of the dynamics of the system, 
a highly simplified model was chosen in [p6| , so 
as not to obscur the relation between the ampli- 
tudes, with complicated interactions and unnec- 
essary quantum numbers. Since this calculation 
does not rely on Eq. (|29|), it can provide addi- 



tional confidence in the correctness of Eq. (32) 
and illustrates how this relation plays out in the 
context of a relativistic field theory. 

Recently, yet another derivation of Eq. (|3^) was 
proposed |2^,^]. The relation between this new 
approach and that of [^6| should be clarified soon 

7.4. Application to the AI = 1/2 rule 

Let n — 1 with W n — ttlk- One possible state- 
ment of the AI = 1/2 rule is that |T°|/|7^| ~ 22. 
What does this rule look like in finite volume? 

For the purpose of illustration, we suppose that 
the scattering phases Sq are accurately given by 
their one-loop expression in xPT JSOj . p| Then, 
using the two-pion energy formulae of Eq. (p9|), 
the size of the box, L, required for the first excited 
state (n = 1) to have energy W\ equal to the kaon 
mass can be calculated in the two isospin chan- 
nels. The results are shown in Table 0. Once the 
size of the boxes is fixed, it is straightforward to 
obtain the proportionality factor which appears 
in Eq. (|32]). The different contributions to this 
factor are also given in Table One ends up 
with 



|T„°| = 44.9x|M°|, 
and 

\T°/TZ\ = 0.92x |M°/M, 

15 Again, when the approach we have laid out is carried 
out in full, these phases will be obtained from a numerical 
study of the two-pion energy spectrum in finite volume. 



I2l=48.7x|^|, (33) 



(34) 



Table 2 

Size of the box, L, required for the first excited, two- 
pion state to have energy equal to uik in the isopin 
and 2 channels. Also given are the contributions to 
the proportionality factor which appears in Eq. (im). 



/ 


L [fm] 


q 


q<t>'{q) 


kSi'(k) 





5.34 


0.89 


4.70 


1.12 


2 


6.09 


1.02 


6.93 


-0.09 



While the factors in Eq. ( j33| ) look large, they are 
mainly due to the relative normalization of free 
states in finite and infinite volume. Indeed, in 
the free case (i.e. with the scattering phases set 
to zero), the relation is \T 7 {\ = 47.7 x |M^|. Thus, 
the effect of interactions is relatively small despite 
the large difference between the scattering phases 
in the two isospin channels (approximatively 45° 
for k = y^m^/4 — m%). This means, in particu- 
lar, that if QCD is to reproduce the AI — 1/2 
enhancement, the large factor will have to come 
from the ratio of finite volume matrix elements, 
\M°/M%\. In fact, as Eq. (fj) suggests, the effect 
should even be slightly larger in finite volume. 

7.5. Summary 

K — > 7T7T rates can, in principle, be obtained 
from the lattice without any model assumptions 
and without analytic continuation pq |. This re- 
quires working on lattices whose sides, L, are 
greater than 5 fm. On these lattices, the effect of 
interactions on the proportionality factor relating 
the transition matrix elements in finite and infi- 
nite volume was found to be small. This means 
that the finite-volume amplitudes must incorpo- 
rate most of the physics which enters the deter- 
mination of their infinite volume counterparts. In 
particular, if QCD is to reproduce the AI =1/2 
enhancement, this enhancement will be clearly 
visible in finite volume. 

For the approach to be fully self-contained, the 
strong scattering phases, which are required to 
relate the finite and infinite volume amplitudes, 
should also be determined on the lattice. As 
shown by Liischer many years ago |57],[52| and as 
briefly discussed above, this can again be done 
using finite-volume techniques. Here, the recent 
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high-precision determination of tttt, S'-wave scat- 
tering lengths by Colangelo et al. [[|| should pro- 
vide a good test of some of the lattice techniques 
required to undertake such studies. 

The same ideas as those described above can 
be applied to baryon decays, such as A — > Ntt, 
X — > Ntt and 5 — > Air, as well as to any other de- 
cay in which the final-state particles scatter only 
elastically. 

There are, however, many potential hurdles in 
implementing the approach. It may be difficult, 
for instance, to extract excited tttt states in prac- 
tice. Furthermore, since the unitarity of the un- 
derlying theory played an essential role in the ar- 
gument leading to the relation of Eq. (|32|), it is 
not clear to what extent this relation holds in 
quenched QCD. Quenched x?T, along the lines 
of p^ , p5[ , should help shed light on this problem. 
Finally, the approach taken here only applies to 
situations where the particles in the final state 
scatter elastically. It says very little about what 
should be done when they can scatter inelasti- 
cally, as for example in B — > tttt, D — > Kn, . . . 
decays. 

8. Conclusions 

This has been an exciting year. New ideas for 
dealing with non-leptonic weak decays and for 
reducing operator mixing with Wilson fermions 
have been proposed. 

There are new chiral perturbation theory re- 
sults relevant for extracting K — > tttt amplitudes 
from K — > 7r and K — > matrix elements and 
from unphysical K — » nir matrix elements. 

There are also new numerical results, many of 
which are still preliminary, for: 

• AS = 2 matrix elements relevant for K°- 

mixing in the standard model and be- 
yond, with domain- wall and O (a)-improved 
Wilson fermions; 

• K — > (tttt)i decays in the CP conserving 
and violating sectors, obtained from K — > 7r 
and K matrix elements, with domain- 
wall (for / = 0, 2) and O (a)-improved Wil- 
son (for 1 = 2) fermions. 

One thus expects, in the near future, lattice 
determinations of Re e'/e and studies of the 



AI = 1/2-rule, hopefully with an active charm. 
In interpreting these results, when they come out, 
it will be important to keep in mind the difficul- 
ties of this calculation and the limitations of the 
K => n and K =r Q approach, as discussed in 
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There should also soon be results for K — > 
(7T7r)/ = 2 decays in the CP conserving and vio- 
lating sectors, from K — > 7T7t matrix elements, 
obtained with O (a)-improved Wilson fermions. 
Kinematical situations other than the traditional 
niK — Tn n with the two pions at rest will be inves- 
tigated. It is further expected that these studies 
will be extended to AI = 1/2 transitions. 

In any event, the coming year may yet be even 
more exciting. 
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